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High-resolution spectral hole burning (SHB) in coherent non-degenerate differential transmission 
spectroscopy discloses spin-trion dynamics in an ensemble of negatively charged quantum dots. In 
the Voigt geometry, stimulated Raman spin coherence gives rise to Stokes and anti-Stokes sidebands 
on top of the trion spectral hole. The prominent feature of an extremely narrow spike at zero 
detuning arises from spin population pulsation dynamics. These SHB features confirm coherent 
electron spin dynamics in charged dots and the linewidths reveal spin spectral diffusion processes. 
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Single electron spin localized in semiconductor quan- 
tum dots (QDs) has attracted a great deal of interest due 
to its potential use in quantum applications il]. Experi- 
mental and theoretical efforts have been focused on con- 
trollable coherent spin dynamics and possible decoher- 
ence mechanisms El H El H [H [TOl El 112 113 El • 
In this paper, we report spectral hole burning (SHB) 
in coherent differential transmission (DT) spectroscopy 
induced by spin-trion dynamics in an ensemble of neg- 
atively charged QDs. Spin coherence induced SHB by 
stimulated Raman excitation and spin relaxation in- 
duced SHB due to the population pulsation dynamics 
are observed in the QD system. Features of SHB not 
only disclose important spin dynamics observed from 
transient spectroscopy [5] and phase modulation tech- 
niques '11] , but also provide information on spin spectral 
diffusion (SD) processes [14, which are not easily revealed 
by previous methods. 

The interface fluctuation GaAs/Alo.sGag.yAs QDs are 
molecular beam epitaxy grown with growth interrupts, 
and modulation Si doping in the barrier incorporates ex- 
cess electrons flE\. The pump Ei{uJi) and probe £^2(^2) 
optical fields are derived from two frequency-stabilized 
and independently tunable CW lasers with a mutual co- 
herence bandwidth of 20 neV, which is crucial for this 
experiment as discussed below. The sample is kept in- 
side a superconducting magnetic liquid helium flow cryo- 
stat and the temperature is maintained at 4.5 K. The 
DT signal is homodyne detected with the probe field by 
a photodiode and extracted by a lock-in amplifier. 

A nonlinear degenerate DT spectrum with uji = loi 
[Fig. [ija)] shows the trion and exciton ensemble reso- 
nances. Their assignments are confirmed by both photo- 
luminescence and transient quantum beats studies (data 
not shown). The trion binding energy (i.e., the separa- 
tion between exciton and trion resonances) is measured 
to be 2.9 meV, in agreement with earlier reports of pho- 
toluminescence [15] and transient spectroscopy ^2 . The 
ensemble trion inhomogeneous broadening width (~ 2.5 
meV) can be estimated from the broad Gaussian profile 




FIG. 1: (a). Coherent degenerate DT spectrum of both 
trion and exciton ensemble resonances. The arrow shows the 
fixed pumping position for (c). (b). Energy level diagram of 
charged QDs in absence of magnetic field, where solid (empty) 
circles indicate electrons (holes) and arrows show the corre- 
sponding spin orientations. The single electron ground state 
has either spin down (|z— )) or spin up (jz-|-)). The trion state 
consists of two-electron singlet state and one heavy hole with 
spin down {\t—)) or spin up (|t-|-)). a~ (cr^) light couples \z—) 
(1 2:+)) to |i— ) (|t+)). (c) A clear view of the trion spectral 
hole burning profile, shown as a double Lorentzian, where the 
probe is detuned relative to pump position (1621.9 meV). 



of the trion resonance. 

To study SHB with non-degenerate DT spectroscopy, 
the pump beam is fixed near the trion ensemble peak 
and the probe beam is detuned (A = a;2 — Wi). A 
narrow spectral structure appears, exhibiting a double- 
Lorentzian-like shape, i.e., a narrower Lorentzian peak 
on top of a broader one [Fig.jljc)]. As established below, 
the narrower Lorentzian peak is due to the trion popu- 
lation relaxation dynamics and its linewidth (^ 12 /ieV) 
gives the trion population relaxation rate. The broader 
Lorentzian peak (^ 50 /^eV) underneath this resonance 
is due to trion population relaxation and coherence de- 
cay broadened by the trion SD process. This complex 
lineshape and its unfolded SHB features with a magnetic 
field in the Voigt geometry are the focus of our study. 

The dynamics of the QD ensemble in the presence of 
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FIG. 2: The calculated SHB lineshape of 2-level system with 
different spectral diffusion rates, where At = lOOOFt and -yt = 
Ft/2, (a), rf^ = 0. (b). rf° = Ft. (c). FP = lOFt. 



the radiation field and various interactions with the envi- 
ronment is described by the modified optical Bloch equa- 
tion (OBE) [IS]: 



dt 



relax 
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where p(e) is the density matrix of each ensemble member 



characterized by its resonance energy e. H is the total 
Hamiltonian including the interaction with the coherent 
optical fields. The second term on the RHS is a general- 
ized relaxation term that describes population decay and 
pure dephasing. The last term is due to various SD pro- 
cesses [T7]. Eq. ([T]) without the SD term reduces to the 
standard OBE. 

Considering the optical selection rules in the absence 
of a magnetic field, the transition from spin ground state 
\z—) (l^-l-)) to trion state \t+) )) is forbidden as a 
dark transition |15j . Therefore, with co-circularly polar- 
ized pump and probe fields, the charged QD energy levels 
reduce to a 2-level system, as shown in the enclosed box 
in Fig.[ljb). The relevant SD process here is the interdot 
transfer of trion population [18] 
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= -rf^(6)af(e) + Wde,e')p,j{e')de', 
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where Wt{e,e') is the spectral redistribution kernel rep- 
resenting the rate for trion population to migrate from 
a dot with resonant energy e' to dots with energy e. 
While this interdot transfer conserves the ensemble trion 
population, the dipole coherence is lost, rf^(e) = 
/ Wt{e' , e)de' is the overall SD rate. For an optically thin 
sample, the DT spectrum is given by the imaginary part 
of the third order induced optical polarization integrated 
over the inhomogeneous distribution |19l I20j , 
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where Fj (jt) is trion population (coherence) decay rate, 
N is the total number of excited charged QDs, Af is 
the ensemble inhomogeneous broadening width of the 
trion resonance and fi the optical dipole of a single dot. 
We have assumed a redistribution kernel Wt{e,e') = 
rf^exp(-(e-e)^/(At)2)/(A/7rAt) where e is the ensemble 
averaged resonance energy [5T]. Eq. ([s]) holds when the 
approximation of plasma dispersion function is taken |22l 
because Ft,7t,Ff^, A <C At and the pump frequency is 
fixed near the center of the ensemble trion spectrum. 

The SD process significantly changes the trion SHB 
lineshape and linewidth, which can be discussed in three 
regimes, as schematically shown in Fig. [2] Ft = 274 is 



assumed only for the theoretical calculation in Fig. [2] 
to simplify discussions, since pure dephasing has been 
found negligible for trion states at 4.5 K [5S]. When 
Tf^ <^ Tt, the trion SHB in the DT spectrum reduces 
to the standard Lorentzian squared with width Fj. 
Secondly, when Ff^ ~ Fj, the trion SHB lineshape is 
(24F2 + A^)[{Tt + Ff + A2]-i[(Ft + 2Ff^)2 + A^]-! 
where the linewidth is considerably broadened by the SD 
process. Finally, when Ff^ ^ Tt, the lineshape changes 
to a double Lorentzian-like profile with larger FWHM of 
2(Ft-|-Ff ^) and smaller FWHM of 2Ft. The physical ori- 
gin of the narrower Lorentzian is due to the population 
pulsation effect [5?^ . This explains the observed lineshape 
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FIG. 3: (a). Energy level scheme of charged dot in Voigt ge- 
ometry {Bx 7^ 0). The new ground states |a;±) denote electron 
spin orientation along the x-axis separated by Zeeman split- 
ting. The new selection rules are labelled by solid (dashed) 
lines with a'^{a~) light, (b). The newly-emerged center sharp 
spike and two symmetric sidebands at — 2.2 T, where the 
fixed pump is indicated by the arrow. The two Stokes and 
anti-Stokes spin-coherence induced sidebands are highlighted. 
The inset shows the zoomed in sharp central spike at A ~ 0. 
A similar figure was first presented in Ref. [11] but with a 
different objective. 



shown in Fig. 



which is consistent with the physi- 



cal model (Eq. (3|) assuming HTt ~ 7fj,eV, H'jt — 8fj,eV, 



46^eV, and HAt 



lOOO^eV plus a slow lin- 
ear slope. The values of Ft and 7t agree with earlier 
reports [21 123] and the value of At is in agreement with 
the degenerate DT trion ensemble resonance in Fig. [ija) . 
The relatively big Tf^ implies trion SD process plays an 
important role. 

With a magnetic field applied perpendicular to QD 
growth direction (i.e., Voigt geometry), the spin eigen- 
states \x±) along the field direction are split in energy 
due to the non-zero electron in-plane g-factor (f/|) while 
the trion states are unaffected as the heavy-hole in-plane 
g-factor is negligible [2S]. The new eigenstates and op- 
tical selection rules are shown in Fig. [sj^a) [2]. With 
co-circularly polarized pump and probe fields, the rel- 
evant states form a three-level A-system enclosed in the 
dashed box. Utilizing a similar A-system, optical pump- 
ing leading to spin cooling has been reported recently in 
different structural QDs [^ ET] . However, no similar op- 
tical pumping effect was observed in our system, where 
the strong spectral diffusion on two spin ground states 
discussed below plays an important role. 

The SHB profile changes dramatically in this Voigt ge- 
ometry. An ultra-sharp central spike appears on top of 
the trion spectral hole, as highlighted by the square re- 
gion in Fig.|3]^b). The physical origin can be understood 
from analysis of the perturbation pathway: 
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The pump field and probe field E2 create a second 

(2) 

order spin population through the population decay 
of the trion, which oscillates at the detuning frequency 



A (known as population pulsations [23]). In the presence 
of the spin population relaxation process, the standard 
OBE predicts a hole burning at A = with a linewidth 
given by the spin relaxation rate ~ F^ [24] . However, as 
shown in Figs. [3| and |4| the measured linewidth is orders 
of magnitude larger than Fg (~ 0.1 neV) measured by 
the phase modulation technique for these QDs [TT]. We 
will establish below that the extra broadening is likely to 
be due to SD processes involving the spin states. 

Similar to the trion SD given in Eq. ([2]), the effects of 
spin SD processes on the spin population and coherence 
are given by. 
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where the density matrix for each QD in the ensemble 
is now characterized by two variables, i.e., the zero field 
trion resonance energy and the spin Zeeman splitting 
in the external magnetic field plus the local field (e.g., nu- 
clear Overhauser field) . Similar to the description of trion 
SD, Ws is the redistribution kernel and Tf^{et,es) — 
J Ws{e't,e'g;et,€s)de[de'g is the spin SD rate. The qual- 
itative feature of the redistribution kernel function de- 
pends critically on the SD mechanism. Local nuclear 
field fluctuation induced SD only affects the spin Zee- 
man splitting: Ws = /(es, ^'s)^{^t — £()■ We note that the 
inhomogeneous broadening of eg induced by the nuclear 
field is given by A^ ~ 0.1 fieV<^ 74 [55]. We also note 
that two quantum dots are equally excited if the differ- 
ence in their resonance frequency for spin to trion transi- 
tion is much smaller than the trion broadening "ft- Thus, 
Px,x{^t,es) = Px,x{e't, 4) if (et - e't) ± (fs - 0/2 < 7t- K 
can then be shown that the two terms on RHS of Eq. ^ 
cancel each other and hence this mechanism has a negli- 
gible effect on the linewidth of the ultra-narrow central 
spike associated with the spin population pulsation dy- 
namics. On the other hand, if the SD process is due to 
the interdot transfer of non-equilibrium spin population, 
it is more reasonable to assume a redistribution kernel 
Ws = rf^exp[-(e, - e-,)VA2 - (ct - e*) VA?]/(^AtA,) 
where et (cs) is the ensemble averaged trion (spin) res- 
onance energy [211 . In the vicinity of zero detuning 
A ^ 7t, the DT signal is determined by 
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shown as the sharp central spike with linewidth of 2rs -|- 
FP, in Yig.^{Theory). The experimental data at vari- 
ous magnetic fields are shown in ¥ig.^{Experiment) , and 
the measured linewidth of the sharp central spike is plot- 
ted in Fig. ^b) from which we extract hVf^ ^ 0.2 /^eV. 
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FIG. 4: Comparison of theoretical calculations and exper- 
imental results in Voigt geometry. Theory- The DT spec- 
trum of 3-level A-system at different magnetic field, where 
At = lOOOFt, 7t ~ rt/2, and rf° = O.Oirt. (a). B ~ so 
that spin inhomogeneous broadening is due to nuclear field 
As ~ A". The central spike and the sidebands merge in the 
limit of vanishing B. (b). B = Bq where fiBdeBo = ShTt- 
Spin inhomogeneous broadening is dominated by the inhomo- 
geneity of the g-factor: As ~ jiB^gtB /h. (c). B — 2Bo. Ex- 
periment: SIfB lineshapes at different magnetic fields, where 
the arrow indicates the fixed pump position (1622 meV) and 
the insets show the zoomed view of central sharp peaks. 



In addition to the sharp central spike, the newly 
emerged SHB features in Voigt geometry also include two 
symmetric sidebands as highlighted by the ellipse regions 
in Fig.jsjb). These sidebands can be understood from the 
perturbation pathway: 
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associated with stimulated Raman spin excitations [2]. 
The Stokes and anti-Stokes sidebands appear respectively 
at A = ±/XB(7gi?/fi, and their separation as a function 
of the magnetic field gives the ensemble averaged elec- 
tron g-factor gl « 0.13 (see Fig. |5ja)), in agreement 
with earlier reports [5]. As the sideband feature is as- 
sociated with the Raman spin coherence, it can be in- 
ferred from Eq. ^ that spin SD process will broaden 
the homogeneous linewidth of the sidebands from 7s to 
7s -I- rf^, where 7^ is the spin decoherence rate. The 
ensemble averaged sideband lineshapes are a convolution 
of the homogeneous lineshape with the inhomogeneous 
broadening of the spin Zeeman energy. In the vicinity of 
A = zLiiBdeB /h, the nonlinear DT signal is 
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At finite magnetic field, the spin inhomogeneous broad- 
ening As(i3) = As(0) -I- HB^deB has the contribution 
from the inhomogeneity of the electron g-factor A^^ in 
addition to the nuclear induced zero field inhomogeneous 
broadening As(0). The sideband linewidth measured at 
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FIG. 5: (a). The separation between Stokes and anti-Stokes 
sidebands is plotted as a function of magnetic field, and the 
linear fitting gives the ensemble averaged electron in-plane 
g-factor \g%\ — 0.13. (b). The linewidth extracted from the 
experimental data at various magnetic fields, where ■ (a) 
are the SHB linewidth of the sidebands (central sharp spike) , 
highlighted in Fig. |3|b) with ellipse (square) region. The solid 
curve is the calculated sidebands linewidth based on Eq. Q 
(see text). The dashed line is a guide to the eye. 



various magnetic fields is shown in FigJSTb), which is con- 
sistent with that calculated using Eq. (pF with parameters 
as h{js + rf^) ~ 0.18 neY and A^ ~ 0.64B /xeV/T -|- 
0.1 /ieV. As 7s <C rf^ from the theoretical investigation 
of both the nuclear and phonon induced spin decoherence 
[HI [13] , we can get nr^D = O.lS^ieV^ which agrees with 
the value of Ff^ extracted from the central narrow spike. 
In addition, the spin g-factor variation Ag^ is determined 
to be 0.01 from As, so Ag^/g^ ~ 0.1, in agreement with 
transient spin quantum beat measurements [2]. 

In summary, we have shown in this paper that SD 
from the trion state complicates the trion SHB profile, 
and makes an important contribution to the double- 
Lorentzian-like lineshape. In the Voigt geometry, we have 
found a complex lineshape arising from spin dynamics in 
the spectral hole burning: a narrow central spike and two 
symmetric Stokes and anti-Stokes sidebands. They have 
been theoretically identified to result as consequence of 
spin population pulsation dynamics and stimulated Ra- 
man spin coherence, respectively. Moreover, a spin SD 
process is observed in contributions to the SHB lineshape, 
and has been theoretically identified as interdot transfer 
of the non-equilibrium spin population. Possible mech- 
anisms include the interdot spin fiip-fiop interactions of 
the electrons or spin conserved electron tunneling to ad- 
jacent neutral dots. Their quantitive estimates can not 
be determined without a detailed calculation which is be- 
yond the scope of the paper. Nevertheless, the existence 
of these decoherence mechanisms revealed by our experi- 
ments will have important impacts on the efforts towards 
spin based quantum applications in these kinds of QDs. 
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